Invariance of the bit error rate 
in the ancilla-assisted homodyne detection 
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We investigated the minimum achievable bit error rate of discriminating binary coherent states 
by homodyne measurements with the help of arbitrary ancilla states coupled via beam splitters. We 
took into account arbitrary pure ancillary states those were updated in real-time by those partial 
measurements and classical feedforward operations where the partial measurements were made by 
homodyne detections with a common fixed phase. It was shown that the minimum bit error rate of 
the system was invariant under these operations. We also discussed how generalize the homodyne 
detection beyond the scheme which made the bit error rate invariant. 
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INTRODUCTION 

Coherent states are fundamental information carriers 
in optical and quantum communications. They transmit 
in pure state intact even through a lossy channel. A 
primitive modulation scheme is the binary phase shift 
keying (BPSK) with | ± a). Because these two signals 
are not orthogonal, i.e. {a\ — a) ^ 0, it has been a 
central issue in quantum communications how to realize 
a quantum receiver that can discriminate them as small 
error as possible. 

Toward approaching the minimum bound for the av- 
erage error probability (the bit error rate, BER), the 
so-called Helstrom bound optimal ^^r9\ and near- 
optimal 0-0 receivers have been studied theoretically. 
Some of them have been put into experimental demon- 
strations [i,[i^. 

These experiments were based on photon counters with 
additional feedback 0, 0, [13] or coherent displacement 
operations However, due to severe limitations on 

their detection efficiency, the results demonstrated so far 
are still far from the theoretical limit without compensat- 
ing the detection efficiencies. Although there is a rapid 
progress in the development of highly efficient photon 
counters [ll|, [l2| , these are still very advanced technolo- 
gies and not always available in quantum optics labs. 

Another important attempt is to use homodyne detec- 
tors. Homodyne detection is already an well matured 
technique, and is used to discriminate the BPSK signals 
in conventional optical communications. Some quantum 
receiver schemes have been proposed in which a homo- 
dyne detector is combined with additional nonlinear pro- 
cesses, such as Kerr medium, or higher nonlinear media 
[3-0] ■ It is, however, still a formidable task to implement 
highly nonlinear unitary processes with low losses. 

Instead of direct nonlinear unitary processes, measure- 
ment induced nonlinear process is another attractive pos- 
sibility to combine with homodyne detectors. Such a 



process is made of non-classical ancillary states, linear 
optical circuit, homodyne detector, and classical feedfor- 
ward system. It was already known that Gaussian oper- 
ations and classical feedforward alone cannot overcome 
the homodyne BER limit |7|]. Therefore non-Gaussian 
elements, or non-Gaussian ancillary inputs, are required 
essentially. Such states are currently available in the lab- 
oratory ISMlTj . It would then be an interesting question 
whether the BER can be reduced below the homodyne 
BER limit by introducing non-Gaussian ancillary inputs 
to the homodyne-feedforward linear optical circuit. 

In this paper, we investigate the achievable BER for 
the BPSK signals when one is available to use arbitrary 
pure ancillary state, beam splitter, homodyne detector, 
and updating the ancillary inputs based on the homo- 
dyne results. We show that under any updating of the 
beam splitter parameters and the ancillary inputs, the 
BER is invariant, being exactly the same as the homo- 
dyne limit. Our result suggests how to set up the system 
with homodyne detection in order to overcome the basic 
homodyne limit. 



MULTI-ANCILLA RECOMBINATIONS 

Let us show our system in FIG.[TJ The received signal 
]s) (s = zba, and without loss of generality a is real) is 
combined with the first ancillary state ]'0i) via a beam 
splitter of a transmittance cos Q\ , the first homodyne de- 
tector, outputting a result y\. This result is used for 
updating the second ancilla and the beam splitter trans- 
mittance. Any explicit expression of the updating rule is 
not needed to show the BER invariance. 

To describe the updating, we introduce notations. An- 
cillary state ]■;/') is expanded by the coherent state as 
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FIG. 1: The received signal |s) is combined with the anciUary 
states \tpi), \tpN) successively via the beam splitters with 
transmittance cos On- Both lipn) and 0„ are updated based on 
the homodyne measurement results of ancillae. The signal is 
finally discriminated by a homodyne detector. The variables 
Xn-i = Un correspond to the quadrature values used for the 
expansion of the transformed signal state in the n-th step. 



The probability density of Eq. ([5]) is 
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where ipicr) = (Q!|V'')j &n.d da = daji A dai for the de- 
composition a = an + iai with aR^ai € R. Denoting 
the signal port (the horizontal line in Fig. [T]) by the sub- 
script "0", and the n-th ancillary port by the subscript 
"n", the whole [n + l)-niode state up to the n-th step is 
recursively defined by 

|*n)o,l:-,n = ■Bo,n(6'n)|*«-l)o,l, ® \'lp,i)n 

d"-a 

—-\s'Jo la'i)iV'i(ai) «"••«) (2) 

The action of the beam splitter Bo,n{d) = exp[0(aja„ — 
aj^fio)] between the 0-th and n-th port is given by the 
matrix representation 
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where s'q = s. Here we notice that the amplitude 
s'^ depends on all the n coherent amplitudes a'-"' = 
(ai, • • • , a„). 

Now let y„ (n = 1, • • • , N) be the n-th homodyne re- 
sult and X be the homodyne result at port "0" . The 
parameter x is a eigenvalue of the quadrature amplitude 
X = (a-|-d^)/%/2, namely (a;|a; = x{x\. The wave function 
is then 

^'Ar(x,y(^)) EE o(a;| i(yi| • • • (g) Ar(yAr|^'Ar)o,... ,N 

d^a 

— • • - VjAf X (a;|s'^)(yi|a'i) • • • {vnWn) (5) 



with 



(x I a) = exp 



(6) 



BIT ERROR RATE WITHOUT FEEDFORWARD 

We first evaluate the BER in the case where the signal 
is discriminated by the {N -t- 1) outcomes of the homo- 
dyne detectors y'''^^ and x without any feedforward in 
Fig. [H Putting xn = x^ we transform the output vari- 
ables {xn,yn) into the new ones (it„,w„) by applying 
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in a descending order n = A^, A^ — I,-- - ,2,1, iteratively 
putting Xn-i = Un- The quadratic parts in Eq. ([7]), say 
is invariant under all of the rotational operations 
(lU, and then is equivalent to 
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Noting that the Jacobian of this transformation is 
1 det i?(6'„) — 1, the conditional probability is now 
expressed in terms of the new variables xo,v*^^) = 
(wi, • • • ,vn) as 

P(a;,yW|s) ^ P(xo,vW|s) = P(xo|s)P(vW (11) 

where 



P{xo\s) = exp -{xo - V2s)'^ 
d^ad^l3 
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Thus it can be factorized into two parts, and the proba- 
bility of obtaining the value xq is independent on those 
ofvW. 
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Therefore the BER can be calculated independent of 
the ancillary homodyne outcomes. So, by introducing 
the threshold 
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determined by P{xth\ct)p{<^) = P{xth\ — oi)p{—a) in a 
priori probability p{±a) , the BER is given by 



dxoP{xo^ v^^'^ \a)p{a) 
+ dxoP(a;o,vW|-a)p(-a) 

•J Xth 



dxoP{xo\a)p{a) + / dxoP{xo\ — a)p{-a) 

= BERo, (14) 

which is identical with the homodyne limit. For the 
equivalent a priori probabilities p{a) — p{—a) — 1/2, 
we have BERq = eiic{V2a)/2 as is well known. 



BIT ERROR RATE WITH FEEDFORWARD 

Next we consider generic feedforward by updating (ipn, 
On) based on the homodyne outcomes y("~i). Suppose 
the parameters (ipmdn) are functions of the previously 
measured values y("~i) of the homodyne detections, and 
then are regarded as functions of wi, • • • , through 
the rotation ([10]). After changing variables via this rota- 
tion (fTO|). the generic functional relations become 



0„(a;o,v("-i)), 



(15) 



for n = 2,3,---,A^. These are shown by eliminating 
y(«-i) from (V'n(y*""^^),6'«(y^""^^)) by iteratively ap- 
plying 

Um = V,n COS 9m + X,n-lSin9m, 

Xjyi X jyi—1 COS Ojyi ^ Til ^lU ^^^^ , 

for m — n — 1, n — 2, • • • ,2,1. 

The parameters {ipi,9i) are given prior to any ancillary 
states. The BERn will be calculated with the help of 
the relation (|15p in the feedforward control case. Note 
that although the probability density of y„ depends on 
the value of signal s, there is no functional relation with 
the parameters {ipn,On) to the value of signal s. The sig- 
nal s affects the value of y„ only through changing its 
probability distribution. The functional relations (|15p 
show that ancillary states do not depend directly on 
the input signal s, but depend only through probabil- 
ity densities of t/„'s. This fact helps us to prove the 



invariance of the BER for the updating operation con- 
sidered here. Because the explicit dependence on the 
parameter s is confined in the probability distribution 
P{xq\s) = exp[— (a;o — V^s)'^]/V^j it is not the variables 
vi,--- ,vn but xq that is affected by the value of the 
signal s = ±a. 

By using the maximum likelihood method, a iV- 
dimensional boundary hyper-surface on the {N + 1)- 
dimensional parameter space {(xo,v(^)) S R^+^} that 
divides the likelihood regions belonging to the values 
s — ±a S R is determined by 

P(xo,v(^V)p(a) = P(a;o,vW| - a)pi-a). (16) 

Here we note that the input signal s = ±a appears in the 
probability P(a;o, v'^^ |s) only through the overall factor 
exp[— (xo — V2s)^], which is the same as Eq. (|12p . Then, 
the boundary condition (fT6|) does not affect the variables 
v(^\ although does determine the threshold point of xo 
as Xq = Xth, which is the same value as that in the sim- 
plest homodyne case. 

Now, let us calculate the bit error rate BER'j^ in the 
feedforward control. The definition 
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is the same as that of BERm except the coefficient ipn 
of the probability density is a function of xq and v^"~^) 
for n = 2, • • • , A^, as seen from Eqs. ([15]). We perform 
the integrations from w^v to Wi in this ordering, using the 
normalization conditions (V'nIV'n) = 1 for n — N, ■ ■ ■ , 1, 
and then integrate xq with the threshold point ([T5)) . 

More concretely, let us consider the first integral J dvN- 
From Eq. ([TS]) the coefficients '02, • • • , V'Af does not de- 
pend on the variable vn, although it does depend on 
the variables wi, . . . , vn-i- Then, we can carry out the 
integration J dwjv by using the normalization condition 
{iPnIiPn) = 1, and we eliminate the coefficient ipN. The 
same step is repeated for f at-i , ■ • • , wi , and then we elim- 
inate the coefficients ■(/'Ar-i,--- ,"01. Finally, there re- 
mains one integral J dxo . Now, this integration becomes 
the same one as that without feedforward and that with- 
out ancillary states (the simplest homodyne case). Fi- 
nally, we obtain the same formula (|14p . 



DISCUSSION 

The invariance proof relies on the transformation prop- 
erty of homodyne measurement base, i.e., the eigenstates 
of the quadrature amplitude and phase. The tensor prod- 
uct of those bases are transformed into a separable state 
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through the beam sphtter 

W)a ® \v')b = Bab{0)\x)a ® \y)b = 

|xcos6' + ysm9)a <8) \ycos9 — xsin9)b, (18) 

where Bab{9) — exp[j0(iaP& — PaXb)]- This is sharply 
contrasted to the beam sphtter transformation property 
of the number state bases. For example, 

Bab\l)a ® |0)6 = C0s6'|l)a ® |0)b + sin61|0)a ® |l)b, (19) 

and the output states are entangled. By the very fact 
of Eq. ([TSl) . the homodyne distributions of the (iV + 1) 
mode systems of the signal and generic ancillae can be 
reduced to the same distribution of the homodyne detec- 
tion without any ancilla. This shows that the homodyne 
detectors fail in causing any useful quantum interference 
to reduce the BER below the homodyne limit. Thus, the 
present setup does not have any necessary non-linearity 
to improve the BER performance. 

The future topic remained is to take into account an 
adaptive updating of the local oscillator phase at each ho- 
modyne measurement. Let us briefly look it at a glance. 
Varying the local oscillator phase in our model is equiva- 
lently realized by generalizing the beam splitting opera- 
tion from the 50(2) group to the U (2) group. One of the 
most general parametrization of the U{2) group elemtnt 
is 
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where <7,j;,ay,az are the Pauli matrices. In the two port 
case of the simplest example in which we have only one 
ancilla and one signal, let us adopt the most general beam 
splitter B{U), the wave function becomes 

1*') = B{Um= I ^\s')®\a'ma), 



U 



(21) 



Next, let us consider the arbitrary phase in the each ho- 
modyne measurement of the signal and the ancilla. The 
phase rotated coordinate base of the homodyne mea- 
surement is given as {x(4>)\q^ = x{x{4>)\ where = 
q cos (p + p sin cf). Then, we obtain 

P{x{cf>o),y{cl>i)\s)^\{x{(t,o),y{Mm\' 

= /^^(a)>(/3)e--(-^")^-(^-"")^ (22) 

where 
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Notice that the phase dependence of a local oscillator, 
say (j), is translated in terms of the phase of the coherent 
state as (x(0)|a) = (a;|ae"^). The rotation of the ho- 
modyne outcomes x, y and the reparametrization of the 
phases of the homodyne measurements forms a group 
50(2) X C/(l)^. If we could reparametrize x,y,(j)o,(j)i 
to eliminate the beam splitter parameter U in the two 
port probability ([22]) , the methodology used in this letter 
would be directly applicable and thus would obtained the 
same conclusion; namely, the BER is invariant. However, 
it is not the case in actual, that is, we cannot absorb all of 
these parameters in Eq. (|20| . Precisely, the unabsorbed 
parameters is in the coset group U{2)/{SO{2) x U{1)^). 
The remained question is to find a way of treating them 
to clarify whether this parameter space helps us to reduce 
the bit error rate further than that of a simple homodyne 
detection. 
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